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The effect of the continuous emission hypothesis on the two-pion Bose- Einstein correiation func- 
tion is discussed and compared with the corresponding results based on the usual freeze-out. 
Sizeable differences in the correlation function appear in these different descriptions of the decou- 
pling process. This means that, when extracting properties of the hot matter formed in high-energy 
heavy-ion collisions from the data, completely different conclusions may be reached according to the 
description of the particle emission process adopted. 
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I. INTRODUCTION 

When describing ultra-relativistic heavy-ion collisions 
with hydrodynamic models, a simple picture has been 
extensively adopted. It is usually considered that, as 
the thermalized matter expands, the system gradually 
cools down and, when the temperature reaches a cer- 
tain freeze-out value Tf, it decouples. Every observed 
quantity is then computed on the hypersurface T = Tf, 
For instance, the momentum distribution of the pro- 
duced hadrons are obtained by using the Cooper-Frye 
integral Q extended over this hypersurface. Though op- 
erationally simple, such a zero-thickness freeze-out hy- 
persurface is clearly a highly idealized concept when ap- 
plied to finite- volume and finite-lifetime systems as those 
formed in high-energy heavy-ion collisions. 

More recently, Grassi, Hama and Kodama Q proposed 
an alternative picture to the particle emission: instead 
of being emitted only when crossing the sharply defined 
freeze-out surface, they considered that the process could 
occur continuously. Being so, in this picture, particles 
could be emitted from the whole expanding volume of 
the system, at different temperatures, and not only from 
the surface with constant T — Tf. As a consequence, 
in the continuous emission model (CEM), the observed 
quantities depend on the whole history of the expand- 
ing system and not only on the instant of the freeze- 
out. Concretely, it has been shown that i) CEM 
enhances the large— my component of the heavy-particle 
(p, A, S, n, ...) rriT spectra, ii) it gives a concave shape for 
the pion mj- spectrum even without considering trans- 
verse expansion of the fluid, iii) it can lead to the correct 
hyperon production ratios and spectrum shapes with con- 
ceptually reasonable choice of parameters ^-Q], and iv) 
it reproduces the observed mass dependence of the slope 
parameter T ||. 

Naturally, we would like to further explore if the above 



model would present striking differences when compared 
to the usual sudden freeze-out picture. One expectation 
would be that the space-time region from which the par- 
ticles were emitted would be quite different in both sce- 
narios. In the continuous emission picture the duration 
of the emission processes is expected to be longer than in 
the freeze-out scenario, which should considerably affect 
the behavior of the correlation function. Previous stud- 
ies have indeed shown that the influence of the emission 
time (6[jj| on the apparent transverse source dimensions 
were remarkably strong. It was also shown in Ref. |||]] 
that a prolonged freeze-out would considerably distort 
the two-particle correlation function. Our main object 
in the present work is to show the differences in two- 
pion correlation predicted by CEM, as compared with 
the results obtained under the usual assumption of sharp 
freeze-out. For this purpose, we will adopt the same ap- 
proximations used in Ref. [|2|, namely, one dimensional 
Bjorken model jIo| for massless-pion gas. It turns out 
that, within these approximations, the HBT effect suffers 
a large deformation when the usual freezeout scenario is 
replaced by CEM, affecting substantially the conclusions 
achieved on the properties of the matter formed in high- 
energy collisions. 



II. CONTINUOUS EMISSION OF PARTICLES 



In CEM, it is assumed that, at each space-time point 
x M , each particle has a certain probability of not collid- 
ing any more, due to the finite dimensions and lifetime 
of the thermalized matter. Then, the distribution func- 
tion f(x, p) of the expanding system has two components, 
one representing the portion of the fluid already free and 
another corresponding to the part still interacting, i.e., 

f(x,p) = ffree(x,p) + fi nt (x,p) ■ 



1 



We may write the portion of free particles as a fraction 
of the total distribution function, as follows 

V 

ffree(x,P) = Vf{x,p) = - _ f in t(x,p) ■ (1) 

Let us assume, as in the previous papers, that the fraction 
still interacting is represented by a thermal distribution 
function 



fint(x,p) « fth{x,p) 



9 



1 



(2tt) 3 exp [p.u{x)/T{x)] ± 1 



(2) 



where is the fluid velocity at x' 1 and T is its temper- 
ature at that point. The factor V can be alternatively 
understood as the probability that a particle with mo- 
mentum p^ escapes from without further collisions. 

If we assume that the fluid is confined to a cylinder of 
radius Rt , the fraction V of free particles at each space- 
time point x^ may be computed by using the Glauber 
formula 



V — exp 



n(x')av re idt' 



(3) 



where 



t out = t + (-pcosc/)+ \J R T ~ /0 2 sin 2 4>)/(vsm9) (4) 

is the time when the particle with velocity v = 
(v sin 9 cos <f> , v sin 8 sin , v cos 9) reaches the surface of 
the fluid at p = Rt ■ 

If we further consider that, initially, the energy den- 
sity is approximately constant (i.e., e = j^T^ for all the 
points with p < Rt and zero for p > Rt), we can calcu- 
late the probability V analytically, resulting in 



V = (t/t 



1.202 



— T Q T(JV re i 



(5) 



where v re i ~ 1. The previous results can be found in Ref. 



III. HBT INTERFEROMETRY 

The second-order interferometry of identical particles, 
also known as HBT effect JllJ is a powerful tool for prob- 
ing geometrical sizes of the space-time zone from which 
they were emitted, as well as for testing dynamical cor- 
relations built in during the system evolution. 

In its idealized version, the two-pion interferometry 
could be studied through the so-called two-particle cor- 
relation function 



C 2 (fci,fe) 



P 1 (k 1 )P 1 {k 2 



(6) 



where Pi(h) and P%{ki,k2) are, respectively, the single- 
particle inclusive distribution and the joint probability 
for detecting two pions; p{k\ — fc 2 ) is the Fourier transform 
of the source space-time distribution. 

In realistic cases, however, it is mandatory to employ 
more general formalisms 12 , as is the case of the 
Covariant Current Ensemble, flexible enough to include 
phase-space correlations resulting from the underlying 
dynamics. As a consequence, the HBT correlation func- 
tions would reflect a model dependent analysis. In the 
Covariant Current Ensemble formalism, the correlation 
function can be expressed as [plJT! 



C(k u k2)=C(q,K) = l + 



\G(q,K)\< 



G(A; 1 ,fc 1 )G(A; 2 ,fc 2 



(7) 



where — kf — k^ and = 2(^1+^2) anc ^ the complex 
amplitude, G(/ci,/c 2 ), can be written as 



G(fc 1; fc 2 ) 



d^x^p e^D^py^hMu^) , 



(8) 



where D(x,p) is the break-up phase-space distribution 
flp|jl2|| and the currents, jo(uf.ki), contain information 
about the production dynamics. If one takes k\ — fc 2 in 
Eq. (||), one obtains 



(9) 



G(k h ki)= / «pfl(i,p)|i„(^ i(1 )i' 



which coincides with the one-particle spectrum. 

As discussed in Ref. the currents jo(uf-ki) in Eqs. 
can be associated to thermal models and written 
covariantly as jo{k) oc y/u^k^ exp {— u M fc M /(2T)}. How- 
ever, to make the computation easier, we shall adopt 
throughout the paper a more convenient parametrization 



u^k 

jo{u.k) = exp{-— ^} 



(10) 



where, in the case of pions, the so-called pseudo temper- 
ature T ps was related with the true temperature T by the 



equation 12 



T ps (x) = 1.42 T(x) - 12.7 MeV 



(11) 



This mapping between T{x) and T ps (x) was later shown 
to be a good approximation also in the case of kaon in- 
terferometry [is}] . 



A. Bjorken model with sudden freezeout 

In the ideal one dimensional Bjorken picture, using the 
above pseudo-thermal parameterization for the currents, 
an analytical form for the amplitudes can be derived E^] 



G(k u k 2 ) = 2<^>{^—J 1 (q T R T )}K a (0 , (12) 
dy q T R T 

where 



S, 2 = [■ypi'm-LT + m 2T ) - ir(miT ~ m 2 T)} 2 + 

2(4^2 +T 2 )miT m 2T [cosh( Ay)- 1] , (13) 

Ay = y\ — y2 and < > indicates average over particles 1 
and 2. 

The single-inclusive distribution is then written as 



^(PN dN ,m iT . 



(14) 



B. Bjorken model with continuous emission 

The initial expectation concerning the differences be- 
tween the continuous emission versus the freeze-out sce- 
narios were mainly focused on the different emission pe- 
riods. Naturally, in the continuous emission picture the 
duration of the emission processes is longer than in the 
freeze-out scenario, which should considerably affect the 
behavior of the correlation function. The reason for this 
comes from previous studies which have shown that the 
influence of the emission time |^-^| on the transverse 
source dimensions were remarkably strong. 

For treating pion interferometry in the case that in- 
terests us, we consider a different but equivalent form 
for expressing the amplitudes in Eq. (Q). The single- 
inclusive distribution is written as in Rcf. 101 



G(ki : ki 



d 4 x [k?f fT 



(15) 



Analogously, the two-particle complex amplitude is 
written, instead of Eq. rt§), as 



G(h,k 2 ) = J d A xe^ x {V^ [k»f free ]}^ {V^ [k%f fr 



]}* • 
(16) 



In Eq.(|15|) and fllq), is the generalized divergence 
operator, which, due to the symmetry of the problem, is 
written in Bjorken+transverse polar coordinates. 

In order to proceed further, let us recall that usually 
we are interested in small momentum differences q 1 * = 
k^ k 2 , as compared with the average momentum of the 
pair, K 11 = ^{k^ + k^). If we then approximate fcf « 
in Eq.(|l6|), a substantial simplification is achieved and it 
could then be written as 



(17) 



G(q, K) ee G(h, k 2 ) = / d 4 x e 1 " x " [K"f f 



We should note that such a dependence on K M , replacing 
the individual momenta k^ and k% in the complex ampli- 
tude of Eq. ( |I6| ) , is also present in the general derivations 
based on the Wigner formalism. 

In principle, the integral in Eq. ( p~7[ ) should be ex- 
tended over the whole space-time with r > tq . However, 
due to the finite size and lifetime of our system, the inte- 
grand is expected to quickly vanish where the assumption 
embodied by Eq. (||) also breaks down. So, in computing 
this integral, we separated the space-time in two regions, 
one where V > and the other with V < , with 
some reasonable value of TV . Upon partial integration, 
the latter is reduced to the surface contribution and the 
former may be estimated by using the Cooper-Frye for- 
mula 0] on the surface V = 7V , applied to the inter- 
acting component. We emphasize, however, that V is a 
momentum-dependent quantity, so this is not the usual 
Cooper-Frye integral. After some manipulation, we get 
for the single-inclusive distribution 



G{ki , k{) 



(2tt)3(1-7V) Jo 



2-7T /> + 0O 



drj 



{ pdpT T m iT cosh(y, - 17) 



r+00 

/ t drpjrkiT cos (j)}e~ 

J Tn 



(18) 



Analogously, instead of Eq. ([L2|), the two-particle com- 
plex amplitude is now written as 



G(q,K) 



1 



(2^(1 -TV) Jo 
t-Rr 



x { / P dp Tjr Mt cosh(y — Tj) 
Jo 



i[T^(qo coshr]—q L sinhr])—pq T cos(0— (j> q )} 



xe 

+00 



To 

x e 



rdrpyrKT cos <j) 

i[r(go coshr;— q L sinh 77) — p^qr cos((p—(p q )] 



X e 



-M T cosh(Y -ri)/T ps {x) 



} 

(19) 



where 

M T = y/ K T + M 2 , K T = + k 2 ) T ,M 2 = K„K" = 

m ' 2 — 3<7m<A Y is the rapidity corresponding to K, <j> is 
the azimuthal angle with respect to the direction of K, 
and 4> q is the angle between the directions of q and K . 

In Eqs. (|l^) and (fjl))), and p? are the limiting val- 
ues corresponding to a certain value of the escape prob- 
ability TV , i.e., 



TF = 



-p cos 4> + v R 2 - ~ P 2 



(kr / E) cosh y sinh 2 (r; — y) + V T 2 ^ a — cosh(r; — y)] 

(20) 



and 



Pt = —t- 



E 



cosh y cos <t>\\] sinh 2 (?7 — y) + V T 2 ^ a — cosh(?7 — y)] 



R*-t 2 < KT ^ 



E 



cosh 2 y sin 2 </> [ y sinh 2 (?7 — y) + V 



-2/a 



sh(?? - y)\ 



1/2 



(21) 



For choosing the value of V? , in principle, we would 
like to take — 1, corresponding to the complete in- 
tegration of ( |l5| ) and (|l7| ). However, we should notice 
that the expressions (|1S|) and ([l9|) above become inde- 
terminate in the limit — > 1. As mentioned above, 
the thermal assumption for fint(x,p) breaks down in the 
same limit. For this reason, already in Ref. g], it was 
chosen V? = 0.5 and the effect of changing this value 
was discussed. We shall adopt the same value here. 
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FIG. 1. 7T7T Correlation Function x qr in Kinematical Zone 
I (side direction), showing the curves corresponding to the 
continuous emission with that corresponding to the sudden 
freezeout. This last case is not sensitive to Kt in this kine- 
matical zone but a slight dependence on Kt can be seen for 
continuous emission. 



IV. COMPARISON OF RESULTS 

A. Ideal Configurations 

The complexity of the expressions for the amplitudes 
appearing in Eq. (j?j) in the continuous-emission scenario 
is evident from Eq.(|I^) and (|l^) above. In order to get 
some insight regarding the differences of the correlation 
functions in the two scenarios under investigation, let us 
select some special kinematical zones, corresponding to 
an idealized situation in which high precision data with 
unlimited statistics would be available. For instance, let 
us fix yi — (Kl = q L = 0), so that 9 = ir/2 with respect 
to the collision axis. Due to the symmetry of the problem 
we can, without any loss of generality, choose K along 
the x-axis. We then explore the behavior of C(qr,KT) 
for fixed Kt- For restricting even more our kinematical 
window, let us consider two cases. Case I (or Zone I) 
corresponds to considering <p pi — —4> P2 — <fi p (the two pi- 
ons are symmetrically emitted around K) , implying that 
4>q = 7r /2; in this case, we can write the individual mo- 
menta as fcf = (y/ml + K T + q T /i , K T , ±q T /2 , 0), 
where the ± signs correspond to pion 1 and 2, respec- 
tively. The momentum difference qr — qs m this situ- 
ation corresponds to the so-called sidewards component 
introduced in Ref For comparison, we consider that 
in the usual freezeout scenario, the decoupling occurs at 
Tf = 170 MeV. The other constant values assumed in 
the calculation that follows were: 
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Results corresponding to the ZONE I above are shown 
in Fig. 1. As expected, since we are neglecting the trans- 
verse expansion, the difference between the predictions 
of the two scenarios is small. Slightly broader correlation 
function for the CEM case, and the decreasing width with 
Kt , was also expected. 

Case II (or Zone II) correponds to considering <p q = 0, 
\ki | > |&2 1 , with both fej along the x-axis, i.e., fcj | 
Kt II qr', in this case, we can write the individual mo- 
menta as fcf = ( yfml + [K T ± qr/2] 2 , K T ±g T /2,0,0), 
where again the ± signs correspond to pion 1 and 2, 
respectively. The momentum difference qr = qo in this 
situation corresponds to the so-called outwards compo- 
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FIG. 2. 7T7T Correlation Function x qr in Kinematical Zone 
II (out direction), showing the curves corresponding to the 
continuous emission with the ones corresponding to the sud- 
den freezeout. 



nent introduced in Ref Results corresponding to 

the ZONE II above are shown in Fig. 2. This is the 
case, mentioned in the Introduction, where the duration 
of the emission process becomes essential Since 
the emission time in the usual freezeout does not de- 
pend crucially on the particle momentum, the correlation 
function is almost independent of Kt ■ On the contrary, 
the emission time is strongly momentum dependent in 
CEM, for large-momentum particles are emitted mainly 
at early times, whereas small-momentum particles may 
be emitted also at later stage of expansion when the fluid 
is cooler and the system larger. So, we see in Fig. 2 a 
significant Kt dependence, being the correlation nar- 
rower for smaller Kt ■ Also, we can see that both the 
CEM curves are narrower than the corresponding freeze- 
out curves, indicating that the emission time in CEM is 
longer in general. Looking more carefully at the curves, 
we can also perceive that the tail of C(qt, Kt) in CEM 
is much flatter than in sharp freezeout. Probably this 
flat tail is due to the small source depth in CEM at early 
times. It is clear that in Bjorken model without trans- 
verse expansion, which we used in the present work, the 
source depth is constant and ~ Rt in sharp freezeout 
scenario. 

Figure 3 represents the same situation as in Zone II 
but with a different freezeout temperature, Tf = 140 
MeV, in order to show the sensitivity of the results for a 
lower Tf in the case of the usual freezeout. As seen, the 
Tf a dependence has shown to be very weak, as expected. 
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FIG. 3. Study of the sensitivity to T/ in the abrupt freeze- 
out model in the Kinematical Zone II (out direction). To- 
gether with previous results shown in Fig. 2, the curve corre- 
sponding to a lower freezeout temperature, Tf a = UOMeV , 
is also included. 



B. Averaged Correlations 

Although the selective kinematical zones could teach 
us interesting points concerning the differences of the be- 



havior of the correlation functions corresponding to both 
scenarios, as shown in Fig. 1-3, such conditions corre- 
spond to an idealization. For putting the calculations 
into more realistic grounds, averages over the angles, mo- 
menta, and the unobserved projections of the momentum 
differences q should be performed. Using the azimuthal 
symmetry of the problem we can still select Kt along 
the a;-axis, such that K = (Kt,0, Kt)- Then, averag- 
ing over different kinematical zones or windows would 
correspond to integrating over K and q (except over the 
plotting component of q). In order to make the analy- 
sis roughly compatible with the range covered by NA35 
S+A collisions [13], we considered the kinematical vari- 
ables in the following intervals: —0.5 < y < 0.5 (or, 
equivalently, -180 < K L < 180 MeV); 50 < K T < 600 
MeV; < {q L lQSiQout) < 30 MeV (corresponding to the 
first experimental bin). As an illustration, we show below 
an example about how to compute the average: 

(C(q L )} = 1 + 

f-Z dK - C dK - Io° Jo 3 ° ^{K, g) |G(*. g)! 2 
JTso dKL IT AKt Jo' d 1 s Jo' d 1oC(K, q)G(k l ,k 1 )G(k2, k 2 ) ' 

The results are presented in the following way. First, 
as done in the preceding subsection, in order to stress the 
differences of results predicted by the two scenarios under 
study we start from the same initial temperature To = 
200 MeV for both the usual freezeout and CEM. The 
results are shown in Figs. 4-6, respectively as functions 
of qL , qo and qs- One sees in Fig. 4 that, as is well 
known, the qL dependence is very sensitive to the freeze- 
out temperature Tf a and if the same initial temperature 
is attained in both scenarios, the correlation function cor- 
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FIG. 4. ( C ) 7T7T x qL averaged over qs and q , showing 
the curve corresponding to the continuous emission hypoth- 
esis with the ones corresponding to the usual freezeout for 
T fo = 170 MeV and T fo = 140 MeV. 



responding to the continuous emission picture is closer to 
the one referring to the thermal freezeout at lower Tf Q . 
However, the shapes are not the same. The one related 
to CEM is more peaked at the small-gz, values, becoming 
flatter in the tail region. This is in clear contrast to those 
corresponding to sharp freezeout scenario which are more 
similar to Gaussians. Physically, this behavior of CEM 
curve could be interpreted as exhibiting the history of the 
matter in expansion, because particles are emitted dur- 
ing the whole evolution in CEM. Namely, the tail of (C) 
depends essentially on early times, when the size of the 
fluid is small and its temperature high, whereas the peak 
reflects later times, when the fluid has fully expanded 
and cooled down. 

We have already seen in the preceding subsection that, 
when plotted as function of qo , the correlation function 
in CEM is significantly narrower than the one in the 
sharp freezeout and has a flatter tail. These features are 
again seen in Fig. 5, where (C) T7r with lower freeze-out 
temperature Tj is closer to the curve for CEM, as in Fig. 
4. However, if the same initial temperature is attainded 
in both scenarios, very low Tf a is necessary, in this case, 
in order to approximately reproduce the same correlation 
predicted by CEM. We can also notice that the depletion 
of the correlation function at small qo values is more dra- 
matic in CEM. In any case, it is important to emphasize 
that our source is totally chaotic in both scenarios and, 
as is well known [||, (C) < 2 at q^ = is originated only 
from the averaging processes, since Coulomb final state 
interactions, as well as the effect of resonances decaying 
into 7r's [^,^| were not considered here. 

We can again notice in Fig. 6 that, also as function of 

Continuous Emission (T ^= 200 MeV) 
Freeze-Out (T = 200 MeV) 
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FIG. 5. ( C ) 7r?r x qo averaged over qL and qs, showing the 
curve corresponding to the continuous emission hypothesis 
as compared to the ones corresponding to usual freezeout at 
Tf = 170 MeV and T fo = 140 MeV, but identical initial 
temperatures To. 
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FIG. 6. { C ) W7T x qs averaged over qL and qo, showing the 
curve corresponding to the continuous emission hypothesis 
as compared to the ones corresponding to usual freezeout at 
T fo = 170 MeV and T fo = 140 MeV, but identical initial 
temperatures To. 

qs , the depletion of ( C ) 7r7r is more pronounced in CEM 
as compared to the curves corresponding to the sharp 
freezeout case. As happened in the previous cases, for 
the same initial temperature, (C) TW with lower Tf a is 
closer to the curve for CEM, but the shape is somewhat 
different. 

In the previous figures, we have shown and discussed 
the differences of CEM correlation function, confronted 
with the usual abrupt freezeout one, when the fluid 
started from the same initial conditions. Now, when ana- 
lyzing the experimental data, the parameters are usually 
adjusted by fitting the data points as close as possible, 
and the conclusions are extracted from the adjusted pa- 
rameters. In the present model calculations, the only 
parameter, besides the freezeout temperature Tf a , is the 
initial temperature Tq. For computing the results pre- 
sented in Figs. 7-9, we have fixed the initial temperature 
for CEM as 200 MeV, and varied the initial tempera- 
ture To for the usual freezeout scenario, trying to get the 
same (or similar) result. To doing so, we have chosen the 
freezeout temperature as Tf Q — 140 MeV, as often done 
in hydrodynamic calculations and following the indica- 
tions of the previous discussions. 

As seen, especially in Fig. 7, the correlation functions 
predicted by the two different scenarios were so different 
in shape that it was not always possible to obtain similar 
curves. For example, in the case of qL dependence, Fig. 
7, (C) for CEM is closer to the curve with higher To at 
small qL, but at high qL, it turns to be closer to the 
one with the lowest value of To. As discussed above in 
connection with Fig. 4, the correlation curve in CEM 
could be interpreted as showing the history of the hot 
matter in expansion, i.e., the tail of (C) reflects essen- 
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FIG. 7. ( C ) lr „ x qL averaged over qs and qo, showing 
the curve corresponding to the continuous emission hypoth- 
esis with the ones corresponding to the usual freeze-out at 
Tf = 140 MeV, but different initial temperatures T . 

tially the early times, when the size of the fluid is small 
and its temperature high, and the peak the later times, 
when the fluid has fully expanded and cooled down. In 
the usual freezeout picture with a fixed Tf a , small To 
is enough to produce a large tail, whereas a larger ex- 
pansion, so higher To , is required to produce a narrower 
(C). 

In Fig. 8, where the qo dependence of (C) wv is shown, 
we can again observe a distortion introduced by CEM 
into the shape of the correlation function. We see that 
(C) for CEM is closer to the curve with T = 230 MeV 
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FIG. 8. ( C )7r7r x qo , averaged over qL and qs, showing 
the curve corresponding to CEM as compared to the ones 
corresponding to the sharp freezeout at Tf a — 140 MeV, but 
different initial temperatures T . 



at small qo, but at high qo, it turns to be closer to the 
one corresponding to T = 260 MeV. However, differently 
from the previous case, Fig. 7, the shape of the freezeout 
correlation curves is only slightly dependent on To and it 
becomes narrower as the temperature increases. What is 
clearly T dependent here is the intercept at the origin, 
which reflects the large expansion dependence of {C) vv 
x q^ as shown in Figs. 4 and 7. 

Finally, Fig. 9 shows that in this case it is possible to 
find an appropriate To for sharp freezeout to reproduce 
the CEM curve. We see that ( C ) for CEM is closer 
to the curve with T = 230 MeV and the agreement is 
good for most of the qs region where the interferomet- 
ric signal is present. This was expected because in the 
present study we neglected the transverse expansion, so 
the transverse size is the same in both the scenarios. The 
initial temperature To in freezeout is higher than the one 
for CEM, because this is required to make the size of the 
fluid large enough and the correlation in the longitudi- 
nal direction sharp enough to decrease the intercept on 
averaging. 



Continuous Emission (T Q = 200 MeV) 
Freeze-Out (T, = 140 MeV) 
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FIG. 9. ( C }nn x qs averaged over qL and qo, showing 
the curve corresponding to CEM as compared to the ones 
corresponding to usual freezeout at Tf a = 140 MeV, but with 
different initial temperature To . 

From the above results, mainly from the correlation 
curves as function of qL , we clearly see deviations from 
the pure Gaussian behavior in cases where the continu- 
ous emission ansatz was assumed. This could actually 
reflect a signature of a continuous process in the particle 
emission. 



V. DISCUSSIONS AND CONCLUDING 
REMARKS 

As mentioned in the Introduction, treating the decou- 
pling process in heavy-ion collisions as occurring on a 
sharply defined surface is an operationally simple but 
highly idealized description. If the consideration of a 
finite thickness of such a decoupling region does not 
bring any noticeable difference in the observable quan- 
tities, such an approximation would be unquestionable. 
However, previous studies |2|-||] have shown that several 
quantities, such as transverse spectra of produced parti- 
cles and heavy-particle production ratios are sensitive to 
more involved description of the process, called continu- 
ous emission model Q. 

In this paper, we concentrated on the two-pion inter- 
ferometry, which has extensively been used as a powerful 
tool for extracting the space-time geometry, as well as 
probing the underlying dynamics of the hadronic matter 
formed in heavy-ion collisions, and studied the differences 
introduced by CEM in confront with the usual sudden 
freezeout. As shown in Sec. IV, also the HBT effect 
suffers a large deformation when the usual freezeout is 
replaced by CEM. This means that conclusions achieved 
on the properties of the matter formed in high-energy 
collisions may differ substantially if we adopt one or the 
other scenario studied here. 

For the sake of conceptual clarity and, evidently, also 
to simplify the computation, we have adopted in this 
work a simplified one-dimensional Bjorken model for 
massless pion fluid (the pion mass has been included only 
to computing observable quantities), without phase tran- 
sition. Nevertheless, one general result emerges, which 
seems to be evident especially by looking at Figs. 7- 
9. Namely, if we describe the same data by using CEM 
or sharp freeze-out (with Tf a ~ 140 MeV), the initial 
temperature To required in CEM is lower than in the 
usual freeze-out. If Tf a is higher, the difference in Tq 
becomes even larger, which is clear from Figs. 4-6. This 
result means that if CEM is the correct description of 
the decoupling process, then it is harder to reach the 
quark-gluon plasma phase than it appears in the usually 
adopted sharp freezeout scenario. 

Since we have worked with a simplified model, we did 
not attempt to make any comparison with data. For do- 
ing this, evidently we have to do some (or all) of the 
following improvements. More realistic equation of state 
(probably including phase transition) should be used; fi- 
nite longitudinal extension of the fluid should be consid- 
ered; transverse expansion should be included; resonance 
formation should be taken into account. All these modifi- 
cations require some hydrodynamic numerical code with 
CEM incorporated. We are now working in this direc- 
tion. 
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